Dipolar collisions of polar molecules in the quantum regime 
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Ultracold polar molecules offer the possibility of exploring quantum gases 
with interparticle interactions that are strong, long-range, and spatially 
anisotropic. This is in stark contrast to the much studied dilute gases of ul- 
tracold atoms, which have isotropic and extremely short-range, or "contact", 
interactions. Furthermore, the large electric dipole moment of polar molecules 
can be tuned with an external electric field; this provides unique opportuni- 
ties such as the control of ultracold chemical reactions a unique platform 
for quantum information processing 121-^], and the realization of novel quantum 



many-body systems 5 



:8|. In spite of intense experimental 



serving the infiuence of dipoles on ultraco 



eflForts aimed at ob- 



d molecules [9|, only recently have 
10|]. Here, we report the experimental 



suflflciently high densities been achieved 
observation of dipolar collisions in an ultracold molecular gas prepared close to 
quantum degeneracy. For modest values of an applied electric field, we observe 
a dramatic increase in the loss rate of fermionic KRb molecules due to ultrcold 
chemical reactions. We find that the loss rate has a steep power-law dependence 
on the induced electric dipole moment, and we show that this dependence can 
be understood with a relatively simple model based on quantum threshold laws 
for scattering of fermionic polar molecules. In addition, we directly observe 
the spatial anisotropy of the dipolar interaction as manifested in measurements 
of the thermodynamics of the dipolar gas. These results demonstrate how the 
long-range dipolar interaction can be used for electric-field control of chemical 
reaction rates in an ultracold polar molecule gas. Furthermore, the large loss 
rates in an applied electric field suggest that creating a long-lived ensemble of 
ultracold polar molecules may require confinement in a two-dimensional trap ge- 
ometry in order to suppress the infiuence of the attractive head-to-tail dipolar 
interactions jlll-ll4l|. 

Dipolar interactions have been explored in several atom gas experiments using the mag- 



netic dipole moments of atoms 



15 



16| . however this interaction is intrinsically orders of 



magnitude weaker than the dipolar interaction between typical polar molecules. Ultracold 
gases of polar molecules, then, open the possibility for realizing strong, and therefore rela- 
tively long-range, interactions. For example, with polar molecules confined in optical lattice 
potentials, one could realize a system where the interactions between particles in neighboring 
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sites is as strong as the on-site interactions now commonly realized with atoms. This longer 
range interaction for polar molecules will allow access to a new regime of strongly corre- 
ate d q uantum gases with unique phase transitions, such as to supersolid phases for bosons 
l3,[l8| and to topological superfiuid phases for fermions [l9|. Another important differ- 



ence between magnetic and electric dipolar interactions is that the strength of the effective 
electric dipole moment is tunable with an applied electric field. In addition to its obvious 
utility for controlling the interaction strength in dipolar quantum gases, the electric-field 
dependence could be exploited in the development of novel quantum computing schemes or 
in the control of ultracold chemical reactions. 

We perform our experiments with an ultracold gas of ^°K^^Rb molecules prepared in a 
single nuclear hyperfine state within the ro-vibronic ground state [lol. ]2o|. The gas 

is confined in a pancake-shaped optical dipole trap, which is realized by overlapping two 
horizontally propagating, elliptically shaped laser beams with a wavelength of A = 1064 nm. 
Typical harmonic trapping frequencies are = 27r x 40 Hz and = 2tt x 280 Hz, in 
the horizontal and vertical directions, respectively. The KRb molecules have a permanent 
electric dipole moment of 0.57 Debye [10], where 1 Debye = 3.336 x 10^^° C m. However, 
the effective molecular dipole moment in the lab frame is zero in the absence of an external 
electric field. When an external electric field is applied, the molecules begin to align with the 
field, and have an induced dipole moment, d, that increases as shown in the inset of Fig. [T]B. 
This figure covers the range of applied electric field that we currently access experimentally, 
which corresponds to an accessible dipole moment range of to 0.22 Debye. In our setup, 
the external electric field points in the vertical direction [z), parallel to the tight axis of 
the optical trap. Thus, the spatially anisotropic dipolar interactions will be predominantly 
repulsive for molecules colliding in the horizontal direction (side-by- side) and predominantly 
attractive for molecules approaching each other along the vertical direction (head-to-tail). 

In an ultracold gas, the quantum statistics of the particles plays an essential role in the 
interactions. Our "^^K^^Rb molecules are fermions prepared in a single internal quantum 
state at a temperature equal to 1.4 times the Fermi temperature. Therefore, the quantum 
statistics requires that the wave function for two colliding molecules be antisymmetric with 
respect to molecule exchange. If one considers the relative angular momentum between two 
colliding molecules, this means that scattering can only proceed via odd partial waves, and 
will be dominated by angular momentum L = 1 (p-wave) scattering at ultralow tempera- 
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FIG. 1: Two-body inelastic loss for fermionic polar molecules. (A) We extract the inelastic loss rate 
coefficient /? from a fit (solid lines) to the measured time evolution of the trapped molecular gas 
density. The data shown here are for an induced dipole moment d = 0.08 Debye (open triangles) 
and d = 0.19 Debye (filled circles), and an initial temperature T = 300 nK. (B) We plot (3, divided 
by T, as a function of the induced dipole moment, d. A fit to a power-law dependence in the 
region d > 0.1 Debye yields a power of 6.1 ± 0.8. The dashed line shows a fit to a simple model 
based on the quantum threshold behavior for tunneling through a dipolar-interaction modified 
p-wave barrier (see text). The solid line shows the result of a more complete quantum scattering 
calculation including an absorbing potential at short range. Note that the full model (dashed line) 
deviates from the simple model prediction at our largest d where the scattering is no longer in the 
Wigner threshold regime. The calculated dependence of the induced dipole moment on the applied 
electric field is shown in the inset to (B). 

ture. Previous work at zero electric field (without long-range dipolar interactions) showed 
that the lifetime of the trapped ^^K^^Rb molecules is limited by atom-exchange chemical 
reactions that proceed via p-wave scattering [2l| . In our experiments, the typical transla- 
tional temperature of the molecular gas is 300 nK, while the energy height of the p-wave 
barrier for '*°K^^Rb molecules corresponds to a temperature of approximately 24 /iK 22 |. 
With the barrier height much larger than typical collision energies, scattering rates in the 
molecular cloud are determined by the tunneling rate through the centrifugal barrier and 
the molecular gas lifetime is relatively long (on the order of 1 s) {21 1. 

In this paper, we investigate the effect of electric dipolar interactions on collisions and 
find a surprisingly large effect even for our relatively modest range of applied electric fields. 
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We measure the molecular loss rate by monitoring the time evolution of the average number 
density of trapped molecules, n (see Fig. [1] A). We fit the data to the solution of 

— = — pn — an, (1) 
dt ' ^ ' 

shown as solid lines in Fig. [1] A. 

The first term on the right hand side accounts for number loss, and we extract the 
measured two- molecule inelastic loss rate coefficient, /3, (which is two times the collisional 
event rate) from the fit. The second term describes density change arising from heating of 
the trapped gas during the measurement. In a single measurement, we observe an increase in 
temperature that is at most 50%. In subsequent analysis, we fit the measured temperature 
to a linear heating and obtain a constant slope c. In Eqn. 1, we then use a = where 
T is the initial temperature of the gas (see also |2l|). 

Fig. [T]B shows a summary of our experimental data, where we plot /3/T, as a function of 
d. We plot the ratio f3/T because the Wigner threshold law for p-wave scattering predicts 
that P scales linearly with T, and we previously verified this temperature dependence at 
d = Debye 2l|. For the data in Fig. [H T ranged from 250 nK to 500 nK. In Fig. [T]B we 
see that dipolar interactions have a dramatic effect on the inelastic collision rate. At low 
electric field, where d is below 0.1 Debye, we observe no significant modification to the zero 
electric- field loss rate (which is plotted at d = 0.01 Debye for inclusion on the logarithmic 
scale). However, for higher electric fields, we observe a rapidly increasing loss rate, with 
well over an order of magnitude increase in /3/T hj d = 0.2 Debye. Fitting the data for 
d > 0.1 Debye, we find that the inelastic rate coefficient follows a power-law dependence on 
d, P/T oc dP, where p = Q.l± 0.8. 

To understand this strong electric-field dependence of the inelastic loss rate, we consider a 
relatively simple quantum tunneling model where the loss is assumed to be due to collisions 
between fermionic molecules that proceed via tunneling through a p-wave centrifugal barrier 



followed by unit probability for loss at short-range [1^. The fact that we do not observe 
any resonant oscillations as a function of E-field (see Fig. 1 B) is consistent with a very high 
loss probability for molecules reaching short range. In an applied electric field, the long- 
range dipole-dipole interaction oc where R is the intermolecular separation, significantly 
modifies the height of the p-wave barrier, and thus changes the inelastic collision rate. 
Moreover, the fact that the dipole-dipole interaction is spatially anisotropic means that 
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FIG. 2: p-wave centrifugal barrier for dipolar collisions between fermionic polar molecules. (A) 
Schematic picture of the effective intermolecular potential for fermionic molecules at zero electric 
field. At intermediate intermolecular separation (center of plot), two colliding molecules are re- 
pelled by a large centrifugal barrier for p-wave collisions. (B) Schematic picture of the effective 
intermolecular potential for a relatively small applied electric field. The dipolar interactions are 
spatially anisotropic, with an attraction that reduces the barrier for head-to-tail collisions and a 
repulsion that increases the barrier for side-by-side collisions. (C) Height of the p-wave barrier as 
a function of dipole moment for repulsive dipole-dipole interactions (Vi, solid red line) and for at- 
tractive dipole-dipole interactions {Vq, dashed blue line). As the induced dipole moment increases, 
the spatially anisotropic dipolar interactions lower the centrifugal barrier for rriL = collisions and 
raise the barrier for rriL = ±1 collisions. The lowering of the tjil = ±1 barrier at very large dipole 
moments is due to mixing with higher partial waves, e.g. L = 3, 5, 7, 

the p-wave barrier height will be different for mi = and rriL = ±1 scattering, where 
rriL is the projection of the relative orbital angular momentum L onto the electric-field 
direction. In particular, the attractive nature of dipole-dipole interaction for polar molecules 
colliding "head-to-tail" lowers the barrier for m/, = collisions, while the repulsive dipole- 
dipole interaction for polar molecules colliding "side- by-side" raises the barrier for rriL = ±1 
collisions. Fig. [2]A,B show these effects schematically, while Fig. |2] C shows the calculated 
maximum height of the mi = and rriL = ±1 collisional barriers, Vq and Vi, respectively, 
as a function of d. 
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In our simple model, we assume that the collision rate follows the Wigner threshold law 
for p-wave inelastic collisions, namely /3 oc T /Vq^^ . For large ci, Vq is significantly smaller 
than Vi, and the loss will proceed predominantly through "head-to-tail" attractive collisions 
of the polar molecules. In this regime, Vq scales as d~'^ and the model predicts that (5/T will 
increase with a characteristic power law of for (i > 0.1 Debye pj:]. This prediction is in 
excellent agreement with our measured dependence of the loss rate on d ioi d > Q.l Debye 
(Fig. 1 B). 

For a quantitative description of the inelastic collisional rate over our full range of ex- 
perimentally accessible dipole moments, we include contributions from both mi, = and 
rriL = ±1 collisions, and we calculate the barrier heights using adiabatic potential curves 
that include mixing with higher partial waves (see Fig. 2 C). We fit the prediction of this 
quantum threshold model to our data using a scaling factor, 7, as a free fit parameter, where 
7 can be interpreted as the loss probability when the collision energy equals the height of 
the barrier. The resulting theoretical prediction (solid line in Fig. [1]) agrees very well with 
our experimental data (open circles), and we obtain 7 = 0.35 ± 0.08. However, in order to 
get this agreement over the full range of d, we found it necessary to introduce a second fit 
parameter that multiplies the coefficient, Cg, of the van-der-Waals interaction. The best fit 
value of this parameter is 1.9 ± 0.9. Also shown in Fig. [T]B as a dashed line is the result of 
a more complete quantum scattering calculation that employs a strong absorptive potential 
at short range but captures the long-range physics and uses Cq as the single fit parameter. 
This fit also agrees well with the experimental data, and gives Cq = 21000 ± 7000 a.u. in 
good agreement with the predicted value of Cq = 16130 a.u. |22[ |. 

Accompanying the increased inelastic loss rates for increasing d, we observe an increased 
heating rate for the polar molecule gas. In Fig. [3] we plot the measured initial heating 
rate as a function of d. The heating rate T = c is extracted using a linear fit to the 
measured molecular cloud temperature vs. time, over a time period sufficiently long to 
allow T to increase by approximately 20 to 30%. In Fig [3], we plot the fractional heating 
rate T/T normalized by both the initial n and the initial T. We have developed a simple 
thermodynamic model for heating that is directly caused by the inelastic loss. We consider 
the energy lost from the gas when molecules are removed in inelastic collisions, and assume 
that the gas stays in thermal equilibrium. In this model (see Supplementary Information), 
the heating arises solely from an "anti-evaporation" mechanism in the trap 23|, |2J] , where 
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Dipole moment (Debye) 

FIG. 3: Normalized fractional heating rate ^/{Tn) as a function of dipole moment. The heating 
rate is extracted using a linear fit to the initial temperature increase and is then normalized by 
the initial density and temperature of the ensemble. The line is the expected heating rate given by 
T/(T^n) = (/3/T)/12 (see text). Typical initial conditions for these data are n = 0.3 • lO^V 

cm 

and T = 0.5 fxK, and the absolute heating rate ranges from 0.1 fiK/s at zero electric field to 2 /iK/s 
at our largest electric fields. 

the particles removed by inelastic collisions have, on average, lower energies than typical 
particles in the gas. One way to understand this heating mechanism is simply to note 
that inelastic collisions preferentially remove particles from the center of the trap where the 
number density is the highest, and where the particles have the lowest potential energy from 
the trap. We also include in our calculations a competing "cooling" effect that comes from 
the fact that the p-wave inelastic collision rate increases linearly with the collision energy. 
Including these two competing effects, we obtain T/(T'^n) = (/3/T)/12 (see Supplementary 
Information). Remarkably, this simple prediction. Fig. |3] (solid line), using the /3/T values 
from the fit to our loss rate data in Fig. [1] with no additional free parameter, agrees very 
well with the independently measured heating rates for the molecule gas. 

The anisotropy of the dipole-dipole interaction is directly revealed in an anisotropic dis- 
tribution of molecules in the trap. The average energy per particle, which we measure from 
the expansion of the gas following a sudden release from the trap, can be different in the 
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FIG. 4: Apparent cross-dimensional rethermalization as a function of dipole moment for > Tx 
(upper row, A-C) and < (lower row, D-F) in the polar molecule gas. The experimental data 
reveal a striking difference between the results for heating the gas in the vertical direction (A to 
C) and heating in the horizontal directions (D to F), and thus provide evidence for the strong 
anisotropic characteristic of dipolar interactions (see text). The electric field is applied along z. 

vertical and horizontal directions. In the following, we present measurements of the time 
evolution of the expansion energy in the two directions for different d. To probe the spa- 
tial anisotropy of dipolar collisions, we start by adding energy along one direction of the 
cylindrically symmetric trap using parametric heating. Here, we modulate the power of 
both optical trapping beams at twice the relevant harmonic trapping frequency for 50 ms 
(z-direction) or 100 ms (x,y-directions). We then wait 100 ms before quickly increasing the 
electric field (in less than 1/is) to the desired final value and measuring the time dependence 
of the vertical (T^) and horizontal (Tj.) "temperatures" of the cloud. Here, and simply 
correspond to the measured expansion energies in the two directions. We note that this type 
of measurement is commonly used in ultracold atom gas experiments to measure the elastic 
collision cross section 25 1. 



Figure H] shows the experimental data for these rethermalization experiments with three 
values of d and under two initial conditions: > (upper row, A-C) and Tz < (lower 
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row, D-F). For d = Debye (Fig. H] A and D), and equilibrate very slowly on a 
timescale of approximately 4 s. Because d = Debye, there are no dipolar interactions and 
the data agree with our expectation of very slow equilibration for spin-polarized fermions. 
Indeed, this is the longest rethermalization time we have observed in our trap, and therefore 
the data are consistent with no elastic collisions and only technical imperfections such as a 
small cross-dimensional coupling in the trapping potential. 

In an applied electric field, the elastic collision cross section due to long-range dipolar 
interactions is predicted to increase as c?^ |26|. For the case where initially > (Fig. 
4 B,C), our data show that and approach each other, in what appears at a casual 
glance to be cross-dimensional rethermalization. The timescale for this apparent rethermal- 
ization even decreases steeply with increasing d as one might expect. However, we note 
that in Fig. 4 C the temperatures cross each other at long times, which is inconsistent with 
rethermalization driven by elastic collisions. Even more striking is the fact that the thermo- 
dynamic behavior of the gas is completely different when the gas initially has < T^. In 
this case, Tz and do not equilibrate during the measurement time (see Fig. |1]E,F). 

The explanation for these surprising observations comes from the spatially anisotropic 
nature of inelastic dipole-dipole collisions and the fact that the molecule gas experiences 
number loss. We have seen (Fig. 3) that loss due to inelastic collisions heats the gas, and 
we can quantitatively understand this heating rate by considering the effect of molecule loss 
on the average energy per particle. We can adapt the heating and inelastic collision model 
described above to allow the average energy per particle, or "temperature", to be different 
in the two trap directions. The "anti-evaporation" picture then predicts that the dominant 
head-to-tail collisions (m^ = 0) will lead to heating in x and y but cooling along z. Side- 
by-side collisions {rriL = ±1), on the other hand, should contribute to heating along z but 
give no temperature change in x and y (see Methods and Supplementary Information). To 
compare this model to our rethermalization-type data, we fix the (i-dependent /3 using the fit 
to our data in Fig. [T] (solid line). This fixes both the time evolution of the molecule number 
as well as the heating rates in the two trap directions. We then accommodate possible elastic 
collision effects in the model by adding a term that would exponentially drive the energy 
difference between the two directions to zero. Figure H] shows a comparison between the 
results of the model (solid lines) and our experimental data. Although the model uses few 
free parameters (only the elastic collision cross-section, in addition to the initial n, T^, and 
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Tz), it provides excellent agreement with the experimental data. 

Our observations provide clear evidence of the anisotropic nature of the dipole-dipole 
interactions through the observed anisotropy in the apparent rethermalization. From the 
agreement of the data with the model, we conclude that the rethermalization behavior is 
actually dominated by the anisotropic nature of the inelastic collisions. What then can we 
conclude about elastic collisions? The best fit value for the elastic collision cross section is 
finite and increases with increasing d. This is consistent with the prediction that the elastic 
collision rate for polar molecules will scale as 26|, and, furthermore, the fit values agree 
with the prediction of cross sections on the order of = 7 ■ 10^'^cm^/(s Debye^) 26|. How- 
ever, the presence of inelastic loss and the resulting anisotropic heating make it difficult to 
accurately extract a measured value of the elastic cross section for a more precise comparison 
with theory. 

The results presented here demonstrate that modest applied electric fields can dramati- 
cally alter the interactions of fermionic polar molecules in the quantum regime. For future 



efforts aimed at studying many-body phenomena in dipolar molecular quantum 
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systems it 
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27|. In 



will be necessary to protect the gas from strong inelastic loss and heating 
particular, the demonstration here of strong spatial anisotropy for inelastic collisions of polar 
molecules suggests that going to a two-dimensional trapped gas will be a promising route to 
realizing a long-lived quantum gas of polar molecules with dipole-dipole interactions. This 
could be achieved, for example, with polar molecules confined in an array of 2D pancakes in 



a ID optical lattice trap 
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12l |. Interestingly, even when short-range inelastic loss processes 



are suppressed, the attractive part of the long-range dipole-dipole interaction could still give 
rise to correlations between neighboring 2D pancakes in the ID lattice trap 28|, |29|. 



A. Methods 

For the fit to a quantum threshold model |l^ in Fig. 1 B, we write the inelastic loss 
rate coefficient, (3 = KqT^ + 2KiT^, as the sum of two terms corresponding to = 
and rriL = ±1 scattering, respectively, and we assume = = T. The (i-dependent 
coefficients, Kq and Ki are obtained using 
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where h = his Planck's constant, and is the Boltzmann constant. The barrier heights, 
Vq and Vi, are taken to be the maximum energy of the long-range adiabatic potential V{R) 
evaluated in a basis set of partial waves \LMl) [14]. The potential V{R) includes a repulsive 
centrifugal term, h'^L{L + l)/(2/ii?^) where /i is the reduced mass of the colliding molecules, 
an attractive isotropic van der Waals interaction —hC^/R^, and the dipolar interaction. We 
use only two fit parameters, h and 7, when fitting this model to the measured /3/T vs. d. 

The solid lines in Fig. 4 are a fit of the measured time evolution of n, T^, and to the 
numerical solution of three differential equations (see Supplementary Information): 

^ = -(i^oT. + 2K,T.^n- - - (3) 

dt V . ^ I X) ^rj.^ \ ) 

^ = ^(-i^oT. + 2iriT,)T, - ^(T, - T,) + c,g (4) 
^ = ^(iroT.)T. + ^-f{T, - T.) + (5) 

Here, we have allowed for a difference in the average energy per particle in the two trap 
directions, "T^" and "Tc", so that (3 = KqTz + 2KiTx. For the fits, we fix the d-dependent 
coefficients Kq and Ki using the previous fit to the inelastic loss rate data in Fig. 1. In 
addition to heating due to inelastic loss, we include a measured background heating rate of 
Chg = 0.01 yuK/s. The elastic collision rate in Eqns. 4 and 5 is given by Tei = ^p^; where the 
elastic collision cross section aei is a fit parameter, v = sj^^^^^^^ ^ind the constant N^ou 
can be thought of as the mean number of collisions per particle required for rethermalization. 
We use Ncoii = 4.1, which was computed for p-wave collisions 30|, however, we note that 
Ncoii depends on the angular dependence of the scattering and may be somewhat different 
for dipolar elastic collisions. 
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Heating due to two-body inelastic p-wave collisions 

Loss of trapped atoms or molecules due to inelastic collisions gives rise to heating through 
a mechanism called "anti-evaporation", where particles removed from the trap have, on 
average, a lower energy than the remaining particles. Here, we model the heating due to 
inelastic two-body p-wave collisional loss in a trapped gas, as relevant to our experiments. 

To accommodate the spatial anisotropy of the dipole-dipole interaction and the possibility 
of an anisotropic energy distribution in the trapped gas, we define "temperatures" along 
three spatial directions using 

where the index i = x,y, z and is the number of molecules. As the trap shape is a fiat 
disk perpendicular to the vertical (z) axis, which is also the direction of the electric field, it 
is natural to assume that the temperatures in the two radial trap directions, Tj, and Ty, are 
equal. The energies Ei are defined as 

^^ = Y1 (^'^^^ + ^^^-^') (2) 

where m is the particle mass, Uz and Ur are the axial and radial trapping frequencies, 
respectively, and the sum is taken over all particles in the trap. The heating rates are then 
given by 

dJ\ ^ _'I\dN IdEi 

dt N dt keN dt ' ^ ^ 

Each two-body inelastic collision removes two particles, so the loss rate is given by 

^ - -^'vr^,. (4) 

where Tcou is the inelastic collision rate per particle. Similarly, 

^ = -AE,NT,oii, (5) 



1 



where AE'j is the average Ei for a pair of particles that is lost from the trap. Putting this 
together, we have 



To find the heating caused by inelastic coUisional loss we need to calculate AEz and AE^. 
The total energy of any pair of particles is the sum of four terms: kinetic and potential 
energy of the center-of-mass motion of the two colliding particles, and their relative kinetic 
and potential energy. In any of the three directions, the ensemble average for the energy is 
2kBT with equal contributions (1/2A';bT) coming from each of the four terms. However, the 
average energy for pairs of particles that undergo inelastic collisions is, in general, different 
from the ensemble average because the collision rate depends on the relative position and 
relative velocity. Because collisions are only sensitive to relative motion, the kinetic and 
potential energy of the center-of-mass motion of the two colliding particles will be on average 
the same as the ensemble average, (l/2fcsT). The relative potential energy is zero for 
colliding particles (since they must be at the same position to collide). This is the reason 
that coUisional loss in a trap leads to "anti-evaporative" heating, in contrast to the elastic- 
collision-based evaporative cooling where energetic particles are preferentially removed from 
a trap. 

It remains then to consider the relative kinetic energy. For p-wave inelastic collisions, the 
collision rate increases linearly with the relative energy of the colliding particles. Further- 
more, as discussed in the main text, the spatial anisotropy of the dipole-dipole interaction 
can result in different rates for ttil = and ttil = ±1 collisions. Therefore, we calculate the 
heating for = and = ±1 collisions separately and then add the heating rates in 
the final result. For uil — collisions, we find the average relative energy for two colliding 
particles using 




coll 



AEi 



(6) 



E\ 



'z,rel 



1 j-oof(Vz,Vr)v'ldv, 

2 JZ^f{v„Vr)vldv, 




(7) 



E' 



'x,rel 



1 S^fjVz, Vr)vl{2'KVr)dVr 
^ !^ f{'^z,Vr){'^T^Vr)dVr 



(8) 



2 



where we assume a Maxwell-Boltzmann distribution of relative velocities given by 

with — vl + Vy. A factor of in both the numerator and denominator of Eqn. (7) 
accounts for the fact that the = p-wave inelastic collision rate scales as collision energy 
oc w^. Similarly, for = ±1 collisions, the collision rate scales as and we find the average 
relative energy for two colliding particles using 

^z,rel = o/" 



2 



2 !^^f{v^,Vr)dv^ 
ksTz (10) 



ksT,. (11) 



Adding this to the ksTi from the center-of-mass energies, we get 

= hnT^ (12) 

A£;° = ^^^t; (13) 

AEl = hsT, (14) 

/^El = 2A;b7;. (15) 

Putting this into the expression for the heating rate (Eqn. (6)), we find that 

^ = -\t^^'coU + \tAoU (16) 

^ = +\t^^Iu- (17) 
In the main text, we defined the loss rate coefficient as 

13 = K^T, + 2K{T, (18) 

where Kq and i^Ti depend on the induced dipole moment. The inelastic collision rate per 
particle is related to (5 through 



_ K^T.n 



co«-^— (19) 

rL« = ^r^- (20) 



With this substitution, our final result for the heating due to inelastic loss is 

"^^^ '':-KoT, + 2K,T,)T, (21) 



dt 4 

^ = ^{KoT,)n. (22) 

If we assume that the ensemble stays in cross-dimensional equilibrium, T — — T^, 
then the overall heating rate is simply 

dT IdT, ^2 dn 



12 

or equivalently. 



dt 3 dt 3 dt 

T) 

{Ko + 2K^)T\ (23) 



^/(T^n) = ^. (24) 
dt'^ ' 12 ^ ' 
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